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ON GENERALIZED KESTEN-MCKAY DISTRIBUTIONS 


PAWEL J. SZABLOWSKI 


Abstract. We examine properties of distributions with the density of the 

form: — r - Ir ^ A 7' y /, C , -r, where c.ai,... ,a n are some parameters and A n 

7rcny =0 ( C ( 1 + a j)- 2a j a: ) 

suitable constant. We find general forms of A n , of k —th moment and of 
k —th polynomial orthogonal with respect to this measure. We indicate their 
connections with polynomials forming the so called Askey-Wilson scheme. On 
the way we prove several identities concerning rational symmetric functions. 
Finally we consider the case of parameters ai,..., a n forming conjugate pairs 
and give multivariate interpretation of obtained distributions at least for the 
cases n = 2,4, 6. 


1. Introduction 


The purpose of this note is to analyze the properties of the following family of 
distributions having densities of the form: 


9 A r n ~ 2 \/r 2 — r 2 

(1.1) fKMKn(x\c, ai ,...,a n )= U V 

7r llj = l( C (l + a U — 2cLjX) 

defined for n > 0, |x| < c, with c > 0, |a ; | < 1. j = 1 Here A n is a 

normalizing constant being the function of parameters ai,,a n . We will call this 
family generalized Kesten-McKay distributions. The name is justified by the fact 
that the distribution with the following density 


( 1 . 2 ) 


fKMK 2 (x\A/a 2 , a, —a) 


v\j4(v — 1) — x 2 
2ir(v 2 — a; 2 ) 


where v = 1+1/a 2 and |a| < 1 has been defined, described and what is more, derived 
in [5] or [9). Then the name Kesten-McKay distribution has been attributed to this 
distribution in the literature that appeared after 1981. 

Thus it is justified to call the distribution defined by ED a generalized Kesten- 
McKay (GKM) distribution. 

Note also that for n = 0 fKMKo{x\c) becomes Wigner or semicircle distribution 
with parameter c. 

It should be underlined that for n = 0,1, 2 distributions of this kind appear not 
only in the context of random matrices, random graphs which is typical application 
of Kesten-Mckay distributions (see e.g.[7], [8], [TO]) but also in the context of the 
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so-called free probability a part of non-commutative probability theory recently 
rapidly developing. One of the first papers where semicircle and related distribution 
appeared in the non-commutative probability context is |fj. 

For n < 5 distribution f kmku can be identified as the special case of the Askey- 
Wilson chain of distribution that make 5 families of polynomials of the so-called 
Askey-Wilson scheme of polynomials, orthogonal. For the reference see e.g. [13] 
For n > 5 the distributions fuMKn were not yet described in detail. 

It has to be noted that in 2009 there appeared paper j3]. Although the aim of 
it was to analyze two dimensional measures on the plane of the form 


V(i-s 2 )(i -v 2 ) 
g(x,y) 


where g was a polynomial in x and y it also contains some results concerning 1— 
dimensional case. The considered there 1— dimensional distribution is very much 
alike the distribution we consider in this paper. The authors of the paper call it 
Bernstein-Szego distribution. The method they use to analyze these distributions 
allow them to consider only the case of even degree n of polynomial g. The results 
of the paper are general and hence rather imprecise. They were obtained by quite 
complicated integration on the complex plane. Our results are precise since we 
assume the precise form of the polynomial in the denominator. Namely we assume 
the knowledge of the roots of this polynomial. Besides as mentioned above the 
order of this polynomial can be also odd. Moreover due to our assumptions we 
are able to give precise form of the polynomials orthogonal with respect to this 
measure. 

Out methods are simple and heavily exploit the properties of Chebyshev polyno¬ 
mials of the second kind allowing to obtain on the way some interesting identities 
like e.g. the ones given in Remark [3] or Corollary [2] 

It should be mentioned that by considering polynomials of even degree, having 
pairwise conjugate roots we are not only able to cover the case of two dimensional 
measures on the plane of the form mentioned above but also generalize the results 
to 3 or more dimensions defying new distributions and finding they marginals. 
Compare Section [2] and the remarks at the end of this section. 

Hence our results give substance and generalize the results of [3j. 

We will present unified approach and recall and collect information on this family 
that is scattered though literature. 

Let us observe first that if X ~ ficMKn{x\c , a ±,..., a n ) and 
Y ~ fKMKn{x 11, a \,..., a n ) then X ~ cY. Hence we will consider further only dis¬ 
tributions with the density fxMKn{x\l, a \,..., a n ) which we will denote for brevity 
fKn{x\ai,.. .,a n ). 

To start analysis let us recall that 


OO 


(1.3) 



where Uk denotes k —th Chebyshev polynomial of the second kind and that 






KESTEN-MCKAY 


3 


Thus the form of the density mu fits the scheme of distributions and orthogonal 
polynomials that was considered in [13; and hence we can use ideas and results 
presented there. First of all let us observe that the densities considered in this paper 
are the special cases of the distributions known from the so-called Askey-Wilson 
scheme of distributions and orthogonal polynomials obtained from the general ones 
by setting q = 0. q is a special parameter called base in the theory of Askey-Wilson 
polynomials. More precisely for n = 0 we deal with Wigner distribution, for n 
= 1 we deal with the so-called continuous big q— Hermite polynomials and the 
distribution that makes them orthogonal. When we set q = 0 and n = 1 we deal 
with distribution fx i (x|ai). For n = 2 we deal with the so-called Al-Salam-Chihara 
polynomials and the distribution that make these polynomials orthogonal. Setting 
now q = 0 leads us to distribution fx 2 {x\ai, 02 ). Further for n = 3 we deal with the 
so-called dual Hahn polynomials and the distribution that make these polynomials 
orthogonal. Finally for n = 4 we deal with the so-called Askey-Wilson polynomials 
and the distribution that make them orthogonal. Hence in particular we know the 
families of orthogonal polynomials that our distributions for n = 0,..., 4 make 
orthogonal. For the precise definitions and further references see e.g. [6] or Cl¬ 
in the sequel we will use exchangeably a„ and (a\,...,a n ) depending on the 
required brevity and clarity. 

Let us denote TI"=i ^ 

We have the following general observation. 

Theorem 1. If aj 7 ^ aj, i ^ j, i,j = 1,... ,n, then 
i) constant A n in m is given by: 


(1.4) 


aT 1 


An — A„(a„) — 1/ y r! 

- %)(1 - aiQj) 


ii) Let us define constants B n ^ for k > 0 by: 


B n , k { a„) = J U k (x)f K n{x\a n )dx, 


Then 


(1.5) 


B n .'k — a u y [ 


^n+k —1 


^ W.fM 1 ' - %■)(! - a i«t) ’ 


for k = 0,.... 

Proof. We will use the fact that 


l ! ~ = 


i -1 


n j= i^(6i - hj) ( x - b i) 
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and the fact that ((1 + a 2 )/a — (1 + b 2 )/b ) = (b — a)( 1 — ab)/(ab). Hence we have 

2A n (-iy 

-V 4 - 

2—1 


* 5 dji) — 


2 n 7r n'li 


- vT^ 2/ - (1 + a?)/(2a*)) 


-\/l —■ 


n—1 TT n 

fl 2 I 1,^2 Oj 


1 


2 A t (-l)* 

2 ^nr=i v " ^ - ot 1 - aw) o - u+a-)/^)) 

2 H n (-l) n /r--A a"' 1 1 


x 1 


-r[ YIj&( a i ~ a i)( l ~ a w) ((! + a v - 2 a i x ) ’ 

Now following (11.31) and orthogonality of polynomials U's with respect to Wigner 
measure we have: 


r 1 2 y / T~-~'E 2 

1—1 tt((1 + a 2 ) - 2<nx) 


Uk{x)dx = a 


□ 


Remark 1. Since coefficients B Ut k are coefficients in the orthogonal expansion 
in ^2 ([—1,1], w), where w denotes measure with the density — \/l — x 2 , in basis 
{Uk}k>o> we 9 e t the following expansion for free: 

r, _ OO 

(1-6) fxnjx |a») = — \/l - x 2 ^ B n ^ k Uk{x). 

7r z ' 

k =0 

For more examples of such expansions see [Cl¬ 


in the sequel we will need the following quantities: 

k 

( 1 - 7 ) S ( k n \a n ) = II a im 

0<ji<j2<-,<jk<nm=l 

(1.8) A,W(a„) = Y, <4 2 ---4-iar 0l+ - +Jn - l) , 

jl + ---+jn-l<m 


for 1 < k < n < m. Whenever it will be obvious we will drop other than k 
arguments of the functions Sk and A*. Notice that S^f' 1 is the k— th elementary 
symmetric function of the variables oi,... ,a n . We set ( a„) = 1. 

Remark 2. We have 


Ax 

A 4 


A 5 

Aq 


3 3 

1 , A 2 = 1 - Oia 2 , A 3 = (1 - aja k ) = (1 - aia 2 )(l - aia 3 )(l - a 2 a 3 ), 

7 — 1 k=j +1 


4 4 

n n (1 - a,ja k )/(l - 54 ) = 
j =1 k=j +1 

(1 - aia 2 )(l - ai« 3 )(l - a 2 a 3 )(l - «i«4)(l ~ Q2Q3)(1 ~ «2a4)(l ~ « 3 a4) 

1 — (21(3.2^3^4 


5 5 

H [] (l- aj a k )/(l-S 4 + SiS 5 -S 2 ), 
j= 1 k=j+l 


rwrwia a i a fc) 

(1 - 5 4 + S 1 S 5 - S 2 -S 6 ~ S 2 Se + S 2 S 6 + S 4 S 6 + S&Se - S 2 - S 2 S 2 + Si)' 
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Notice that for n = 1,... ,4 the A n agrees with respective constants given in |l3j 
in formulae (2.f), (2.6) and unnamed formulae on top of 9-th and 10-th pages for 
q = 0, when presenting respectively densities that make big q-Hermite, Al-Salam- 
Chihara, dual Hahn and Askey- Wilson polynomials orthogonal. 

(2) (3) (3) 

Using this denotations we have: B 2 ,k = A), , B 3 ^ = A),' — S 3 Compare 
here with formulae (2.7) and (2.10) in [13] for q = 0. 

Remark 3. Let us notice also that 

/ I n _ n 00 

Ai/l - x 2 JJ E a l f i Uk i (x)dx = 

_1 i=l ki=l 

n 

ki> 0 ,...,k n >0 j =1 

where 14 1 ,...,k n = f\ f t/1 — x 2 nT=i Uki(x)dx. That is we get generating function 
of the numbers Vk lt ... t k n for free. This observation gives also interesting interpreta¬ 
tion of constants A n . 


Corollary 1. i) 

n 1 ^ L^J /^ 

J ^x k f K n(x\a n )dx = + 1 ^ 9fc E ( fc - 2 J + 1 \ j J B n,k-2j{a„), 


a) 


min(m,A;) 

Uk(x)U m (x) fKn(x\a n )dx = E -Bn, \m—k\+2j («*n) • 
1 i=0 


Proof, i) we used well known identity: 2 k x k = — (j-i))Uk- 2 j(x) (see 

e.g. [S] Proposition 1 with q = 0 and the fact that h n (x |0) = U n (x )) and the fact 
that ( k ) - (E) = (k- 2 j + 1 )( k+ j 1 )/(k + 1). ii) we used identity U k (x)U m (x ) = 
^mm(Avm) (a;) that can be easily derived from 2T n {x)T m (x) = T n+m (x) + 

T\ n -m\ (x), where T n denotes Chebyshev polynomial of the first kind and the formu¬ 
lae relating Chebyshev polynomials of the first and second kind. See also formula 
(2.13) of [E] with q = 0. □ 


Proposition 1. For any function g : R — > R let us denote g(a n ) = (g(a 1 ), ..., g(a n )). 
Also let = ( 61 ,..., bi-i, 4+ 1 ,..., b n ). Let g(x) = (l+x 2 )/(2a;). Then if ai 7 ^ aj, 
i y£ j, i,j = 1 ,... ,n, we have for every k = 1 ,... , n — 1 : 


(1.9) 


y n <~ 2 Sk( gtQ) 

E' n j#,( a i - a i)(! - a i a i) 


= 0 . 


Proof. We start from the fact that 1 / — 4) = 5Z”=i D n i 7 w ^ ere -Dn j 
= viw — bj). Thus from the properties of simple fractions decompositions 
we have the following identity 

n n n 

!/ HE ~ = (E Dn IK* ~ IK* _ 

i —1 *=1 i=l 
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Hence all coefficients by nonzero powers of x must be zero. In particular for Vfc = 
0 ..., n — I we must have Y17=i = 0. Now it remains to substitute 

= g(a,i) and use the fact that (&,; — bj) = ( aj — ctj)( 1 — ciiCij)/( 2 aiCij) so that 

n n 

A»,» = 2 n_1 a"" 8 (JJ a 0 )/ JJ((1 - - a,)) 

i= 1 


lj5« a 3 ~ a i Ely = 1 a j- 


smc ea- * flj 

As a corollary we get in particular the following identities: 

Corollary 2. Vn > 2 : i) 


□ 


n-1 


(f.10) 

ii) 

Hi) Vto > 1, n > 2 : 
( 1 . 11 ) 


y 

^ n ^i( a i - aj)( 1 - «i«i) 

" a r 2 (i+a?)E J#< n fe# ia fc 


= 0, 


= 0, 


an) = 0. 

i=0 


Proof. To prove i) we take fc = 0 in m- To prove ii) we take k = 1 in (11.91) and 
notice that 


n— 1 


(II a i) S 'i(s( a ! 


(<h = <-^(1 + ^ 11 ^ = 

n—2 ( 


«r 2 £Il“‘+«"En 

iV» Mi iV* Mi 

To prove iii) we observe first that n > 1, fc > 0, Vi = 1,... n : 


a k = af '(l + a 2 )y ]Ja fc . 

M* Mi 


( 1 . 12 ) 


Pa 7 +k - j Sf{ a n ) = 0, 


i=o 

which is elementary to notice. Secondly, using «E 3 we have: 

Tt 71 71 T7, - J - 7T1 /u 1 

y (-i) fc ^ n) H n „ m _ fc = A n y (-i) fe ^ n) y ■ 


/c—0 


k—0 


n - aj)( 1 - aittj) 


1 


- • 4 ” S n;t>. 4x1-■,«,) - 0. 


by (11.121) . 


□ 


Theorem 2. For every n > 0 i/ie family of polynomials orthogonal with respect to 
f Kn is of the form: 

nA(2ra+2) 

Pi") (,r|a„) = y (—1)1 U m -j (x) Sj 1 ^ (a„), 

i=o 


(1.13) 








KESTEN-MCKAY 


7 


to > 0, where is given by 

Proof. Since order of the polynomial P m must be equal to to and we have the 
following identity : 

U—m—2i x ') — blm{ x \ 

both those facts explain why we must have A (2to + 2) as upper summation limit 

in (irm 

(n) 

The fact that polynomials Pm can be expressed as linear combination of the 
last at most n + 1 polynomials U m follows directly from H] (Proposition 1 iii)) or 
from [3] (Lemma 3.1). Similar fact was noticed for n = 1,2 earlier by Maroni in 
his papers published in the 90’s. However the form of coefficients of this linear 
combination can be most easily found arguing by induction. 

Assertion of the theorem is true for n = 0, 1. For the proof see e.g. formula 
(14.18.4) of [6] with q = 0 or e.g. pE5] (3.1 with q = 0). We see that then h m (x\a, 0) 
= U m {x) — aU m -i(x), to > 1 and h 0 (x\a, 0) = 1. Now recall Theorem 1 of pREj. It 
states that if we have two measures dA and dB related by the formula dA = -Af^dB 
and we know polynomials orthogonal {b m } with respect to dB then there exist a 
numerical sequence n m such that to— th polynomial d m orthogonal with respect to 
dA is equal to d m = b m + K m b m _i. By the way the sequence K m is defined by the 
recursion:«; m = /3 m _ 1 — /3 m _i/K m -i + D , with n\ = /3 0 + D — C where sequences 
(3 m and $ m are the coefficients of the 3-term recurrence satisfied by polynomials 
b m . More precisely b m+ \(x) = (x — /3 m )b m (x) — $ m _ 1 b jn .~i(x). Hence by induction 
assumption we have 

nA(2m+2) 

p£\x\a n ) = ( — 1 y Um-j (x)Sj H ' > (a„ ) . 

1=0 

Further by the quoted Theorem we know that 

(1.14) P^n +1 \x |a„+i) = P^\x |a„) + k^ + 1) P^l^x |a„). 

Notice that from (fL6l) and (11.111) it follows that to > 1 : 

J Pm\ x \ a n)f K ( n )(x\a n )dx = 0 . 

Now let us notice that constant Km +1 ^ can be found from the condition 
J Pm +1 \ X \ a n+l) fK(n+l)( x \ a n+l)dx = 0 

for all to > 1. 

Using (11.141) we deduce that 

raA(2m+2) nKlm 

p£ +l \ x |a n+1 )= £ (-lyUm^ws^M + K^+v £(-i 

1=0 1=0 

= U m (x) - U m -i(x)(s[ n ^ (a„) - k<£ +1) ) + U m _ 2 (x)(S^ n \a n ) - K f 1 )s{" ) (a n ))... 

Now we see that for k|I +1 * = — a n+ i we get S[ n \sLn ) — k„ + 1 ' = S'J” +1 ' ) (a ra+1 ), 
Seisin) - Km +1 ' > s[ n \ a n ) = S , 2 ” +1) (a n+ i) and so on. Finally Km +1) Sn\ a„) = 
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—5 , ^^ 1 ( a „ + i). Hence with such choice of Km +1 * we get the following formula 


(n+1) A(2ra+2) 




1. □ 


Remark 4. Notice that polynomials Pm\x |a„) are neither monic nor orthonor¬ 
mal. 


Remark 5. Notice that for m > n we have the following 3 term recurrence formula: 


2xP m (x| a r/ ) — Pm+i(x|a n ) -+■ P rn —i(^x |a n ), 

which follows directly similar formula to be satisfied by the Chebyshev polynomials 
{U n } . For m < n the problem is more complex. It follows the fact that we naturally 
have:P-i(x\a n ) = 0, Po(x|a„) = 1. However we also have 

Pi(x|a„) = 2x(l - 5 4 (a„) - 5i(a„) + 5 3 (a„), 


which results in 

Pi(x|ai) = 2x — ai, Pi(x|a 2 ) = 2x — (oi + a 2 ), 

Pi(x|a 3 ) = 2x — <Si(a 3 ) + S^af), Pi(x|a 4 ) = 2x(l - S , 4 (a 4 )) - Si(a 4 ) + S 3 (a 4 ), 

and so on. Note that for n = 1,2,3,4 we get polynomials of order 1 of Askey- 
Wilson scheme for q = 0, i.e. big-Hermite (n = 1), Al-Salam-Chihara ( n = 2), 
continuous dual Hahn (n = 3), and Askey-Wilson (n = A). 

Similarly for m = 2 we have 
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P 2 (x|a„)=^(-l)^ 2 _,(x)S> n ) 


which gives 

P 2 (x|ai) = U 2 (x) - aiPi(x), P 2 (x|a 2 ) = P 2 (x) - (a 4 + a 2 )U 1 (x) + aia 2 , 

-P 2 (x|a 3 ) = t/ 2 (x) - Pi(a 3 )Pi(x) + ^(as), P 2 (x|a 4 ) = P 2 (x) - S'i(a 4 )Pi(x) + 5' 2 (a 4 ) - ,S' 4 (a 4 ), 

Similarly we notice connection with big-Hermite (n = 1 ), Al-Salam-Chihara ( n = 

2), continuous dual Hahn (n = 5), and Askey-Wilson (n = 4) polynomials for 
q = 0. Compare [T7] (2.3 with set q = 0 and successively a = d, d = O&c = 0 ,d = 

O&c = 0&6 = 0), 

Proposition 2. Let X ~ fKn{x\a n ) then V |t| < 1 : 


j-p _ '*VTl\ v \ TL) 

1 + t 2 - 2 tx = n?=i(l-tai)’ 


1 _ Qn(t |a n ) 


where 



is a polynomial of degree max(n — 2, 0) in t. 
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Proof. Recall that Y^’jLo t n U n (x) = 1/(1 + t 2 — 2tx). Thus 


r*i °° 


00 n 


a. 


n+j-l 


(y Unjx)) fKn{x\&7i)dx — A n ^ ' P ^ ' y-j-n , w-, \ 

]^o n&i(Oi - Oj) (! - a ^ a o ) 


_1 J=0 


= ^«E 


,n-1 


jrf 11,- a i)( 1 - a i a j) 


E*' a * = A»Et~ 


n—1 


i>o 


^ 1 - a it n - «i a j) 


nr=i(i-tai 


;AE»r'n 


(1 - a,t) 




(a* - aj)(l ~ a i°i) 


The fact that Qn is a polynomial of degree n — 2 follows the fact that a/ 

CLjt) is a polynomial of degree n— 1, but the coefficient by t n 1 is equal to S n ( a„)a" 2 ' 
The assertion follows from (11.101) . □ 


Remark 6. We have by direct calculations: 

Qi(t\a) = l,Q 2 (t\a 2 ) = l,Q 3 (t|a 3 ) = 1 -tS 3 (a 3 ), 

Qi(t |a 4 ) = ((1 - S 4 ) - t(S 3 - S 1 S 4 ) + t 2 S 4 ( 1 - S 4 ))/(l - S 4 ). 

As an immediate consequence of this formula and definition of B n ^ we get 
characteristic functions of numbers -B n ,fc(a„): 

Corollary 3. Vn > 0 : 

n 

E tkB n,k( a n) = Qn(t\a n )/ J|(l - taf). 
k> 0 i= 1 


2. Complex parameters 

In this section we will study properties of the generalized Kesten-MacKay dis¬ 
tributions for even n = 2k and parameters ai , i = 1 ,... ,2k being complex and 
forming conjugate pairs. The new parameters will have new names namely for the 
conjugate the pair say a* = p. t exp [iOf) and = p t exp(— iOf) we will denote yi 

= cos 9i for i = 1,..., k. Besides notice that we have: 

(1 + p 2 exp(2 i6i) - 2xp i exp(*0,))(l + p 2 exp(-2 iSf) - 2xp i exp(—*0*)) 

= 1 + pf + 4:X 2 p 2 - 4xp,(l + p 2 ) cos 9i + 2p 2 cos 29 i = w(x,y i \p i ), 

where we denoted for simplicity: 

(2.1) w(x, y\p) = (1 - p 2 ) 2 - 4:ryp(l + p 2 ) + 4p 2 (a; 2 + y 2 ). 

Hence now the density fK 2 k(x\a\,..., a 2k ) will have the following form that we 
will denote by fMk(x\y k , p k ) : 


( 2 - 2 ) fMk{x\y k ,p k ) = Mk — r -;-—, 

71 11 J = 1 W \ X 1 Ui\Pi) 

with I Pi\ < 1 , \yf\ < 1 . 

Following Remark [2] we have: 

Lemma 1. i) A 2 = 1 — p 2 , 

ii) S[ A \a 4 ) = yipi+y- 2 p 2 , -^(a*) = p\+ pl+^ymPxP^ si 4 \a4) = 2p l p 2 (p 1 y 1 
+ p 2 y 2 ), A4 = (1 - pf)(l - P 2 )w(y 1 ,y 2 \p 1 p 2 )/(1 - pfpl) 
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in) s[ 6 \&q) = p x yi+ p 2 y 2 + P 3 y 3 , S^\a 6 ) = P? + p| + p1 + 4(piP 2 2/i?/2 + P 1 P 3 P 1 Z /3 

+ P2P3V2V3), 

4 6) ( a 6) = 2 (p? + pl)p 2 V 2 + 2(p| + pDp-lV 1 + 2(p 2 + p 2 2 )p 3 y 3 + + 8 p 1 p 2 p 3 y 1 y 2 y 3 , 

Si 6) { a 6 ) = PiP| + p|p| + Pip! + 4p 1 p 2 p 3 (p 3 yi2/ 2 + P 2 yiy 3 + Pi 2 / 22 / 3 ), 

= 2 PiP 2 P 3 (PiP 2 V3 + PiP 3 V 2 + P 2 P 3 yi), S^ e \ a 6 ) = p?plp§, 

A 6 = (l-p0(l-p|)(l-P3) w (Pl ) y2|PiP 2 )w(p 2 ,?/3|P2P3) u; (2/l I 2/3|PlP3)/^3(?/i,?/2,2/3,Pl,P2,P3)» 
where we denoted 

(2.3) w3(y 1 ,y 2 ,y 3 ,p 11 p 2 ,p 3 ) = (1 - p?p|)(l - P 2 p1)(1 ~ P 1 P 3 X 1 - P 1 P 2 P 3 ) 

-4piP 2 P 3 (l + P 1 P 2 P 3 )(Pil 1 - P 2 X 1 - pl)y2V3 + P 2 ( 1 - P?)(l - P 3 )2/i2/3 + P 3 (l - P?)( 1 - P 2 ) 2 /l 2 / 2 ) 
+4p?p|pi((l - p?)( 1 - P 2 P 3 )yl + (1 - pl)(l - PiP 3 )2/i + (i - P§)(1 - P?P^)P 3 2 )- 
Proof. All calculations were done using Mathematica 10. □ 


Remark 7. Recall that for all |x| , \y\ < 1 and \p\ < 1 we have useful expansion: 

1-p 2 


(2.4) 


'Ax,y\p) 


= Y / P J U J (x)U J (y), 
i=o 


which is nothing else but the famous Poisson-Mehler formula for q = 0 (for the 
reference see e.g. 0) (13.1.24) or for alternative proof [11 ] ) 


Following the above mentioned remark we have: 

(2.5) 

9 a _ 00 k 

fMk (x I y fc ,p k )= \/l-x 2 ^ 11 /'V' 1 (•< )/ , fij (// y ). 

Pj) mi 0 j=l 

k 00 k 

(2.6) ]J(1- p 2 )/A 2k = V mi ,..., mk '[[p?U mj (y j ), 

j =1 mi,m2,...mfe=0 j =1 

where as before, above 14 1 ,...,k„ = f _ 3 — x 2 Y\a=i Uk i {x)dx. 

Since each density fMk can be presented as linear combination of fM\{x\yi, pf) , 
k = 1,..., k (by simple fraction decomposition) we will analyze /mi first. We have 
the following result: 

Theorem 3. i) My G [-1,1] : f\ f M i(x\y, p) = 1, 
ii) f\ fMi(x\y,p)ly/l - y 2 dy = f\/l - x 2 , 

Hi) f_ 1 fMi(x\yi,p 1 )f M i(yi\y 2 ,p 2 )dyi = f M i(x\y 2l p 1 p 2 ). 


Proof, i) Either we use directly properties of fK 2 {x\a,b) and the fact that in our 
case ab = p 2 , or we apply (12.41) . ii) follows directly from i) and the fact that 
w(x,y\p) = w(y,x\p). iii) we have: 

/ /mi(£|2/i,Pi)/mi(2/i|2/2,P 2 )^2/i = 


2 
7T 
2 
7r 


Vl~x 2 J 


1 2 {i - pD(i ~ pDV 1 - y‘i 


_! 7T w(y 1 ,x\p 1 )w(y 1 ,y 2 \p 2 ) dVl 
-(! - Pi)(l - pl)Vl - x 2 /A 4 = f M 2 {x\y 2 ,p 1 p 2 ), 


since A 4 = (1 - p?)(l - p%)w(x, y 2 |PiP 2 )/(l - pIpD- 


□ 
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Remark 8. Results of the above Theorem indicate possible applications of the dis¬ 
tributions /mi and Wigner in multivariate analysis and stochastic processes. More 
precisely assertion i) shows that fMi(x\y, p) is in fact a conditional distribution, 
ii) shows that fM\{x\y, p) f mo(v) can be treated as a density of certain bivariate 
distribution with /mo that is Wigner distribution as its marginals, finally Hi) is 
nothing else but the so-called Chapman-Kolmogorov property. These properties are 
known and applied in stochastic processes, see e.g. [2] and |12] . We quoted them 
for the sake of completeness of the paper and also in order to present new proofs 
of these properties directly basing on the general properties of generalized Kesten 
distributions discussed in the first part of this paper. 


Following the above mentioned remark let us denote by f 2 (x,y\p) denote 2— 
dimensional measure defined by: 


(2.7) f 2 M(x,y\p) = fMi(x\y,p)f M0 (y) = 
Remark 9. Notice also that 

fM2{x\yi,p 1 ,y2,p 2 ) = 


(1- P 2 )V(l-x 2 )(l-y 2 ) 
ir 2 w(x, y\p) 

2 Vl - x 2 (l - p\){ 1 - pl)w(yi,y 2 \p 1 p 2 ) 


7T w(x,y 1 \p 1 )w(x,y 2 \p 2 ){l - pjpl) 
fMi(yi\x,p 1 )fMi(x\y2,P2)fMo(y2) 
fMi(yi\y2,PiP 2 )fMo(y2) 

which can be interpreted in the following way. Let us consider 3 element discrete 
Markov chain X 3 ,X 2 ,X 3 such that transition density X 2 \X 3 is fMi(x\y 2 , p 2 ), tran¬ 
sition Xi\X 2 is fMi{x\y\, pi) while marginal density of X 3 is /mo ( 2 / 2 ) then the 
conditional density of X 2 \X 3 ,X 3 is fM 2 {x\yi,Pi,y 2 ,p 2 ). 

w 3 (yi, V2, 1/3, P1.P2.P3) 


Lemma 2. f\ f-^1 - y\ 


w(yi,y2\p 1 p 2 )w(y 2 ,y3\p 2 p 3 )w(yi,y3\piP3 ) 


dy 1 = 


1 2 2 

1~P 2 P3 


w(y2,y3,\p 2 p 3 ) ‘ 


Proof. We start from (E3D considered for k = 3 and get n;=i(i - pD/m = 


_ w3(yi,y 2 ,y3,Pi ,p 2 ,p 2 ) _ 

■u’{yi,y 2 \p 1 P 2 )w{y 2 ,y 3 \p 2 p 3 )w(y 1 ,y 3 \p 1 p 3 ) 

E oo v^oo n mi n rn 2 n rn3 V 

m i=0 Z^m 2 = 0 Z^m 3 =0 Pi P 2 P 3 v r. 


mi 77127723 


Umi {Vi ) u m2 {y 2 ) u m3 (y 3 ) basing on Lemma 


ED Now f 1 1 (yi.y^a.P^.P,) -- d = 

I—I J- 7T V *1 w{y 1 ,y 2 \p 1 p 2 )w{y 2 ,y 3 \p 2 p 3 )w{y 1 ,y 3 \p 1 p 3 ) y 

Em 2 =oEm 3 =oft 2 P's 3 ^0m 2 m 3 tlm 2 (j/2) tlm 3 (?/3) = EEm 2 =0 {P2Ps) ^ Um 2 (j/2) U m2 (?/3) 


P 2 P 3 


□ 


w{y 2 ,V3.\p 2 p 3 )' 

From this Lemma we derive the following important conclusion. Namely that 
the function: 


( 2 . 8 ) 


g{yi,y2,ys\pi,P2iPs) = 




_ w 3 (y 1 ,y 2 ,y 3 ,p 1 ,p 2 ,p 3 ) _ 

w(yi,y 2 \p 1 p 2 )w(y 2 ,y 3 \p 2 p 3 )w(yi,y 3 \p 1 p 3 ) 


can be treated as the density of some 3 D distribution with 2D marginals equal to 

f2M(yi,V2\pip2), f 2 M(yi,y 3 \PiP 3 ), /2m(j/2, ys\p 2 p 3 )- 
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